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By performing first principles calculations we investigate the intrinsic anomalous Hall conductivity 
(AHC) and its anisotropy in ordered Llo FePt, CoPt and NiPt ferromagnets, and their intermediate 
alloys. We demonstrate that the AHC in this family of compounds depends strongly on the direction 
of the magnetization M in the crystal. We predict that such pronounced orientational dependence 
in combination with the general decreasing trend of the AHC when going from FePt to NiPt leads 
to a sign change of the AHC upon rotating the magnetization direction in the crystal of CoPt alloy. 
We also suggest that for a range of concentration x in CoxNii-^Pt and Fe^Coi-ajPt alloys it is 
possible to achieve a complete quenching of the anomalous Hall current for a certain direction of 
the magnetization in the crystal. By analyzing the spin-resolved AHC in 3dPt alloys we endeavor to 
relate the overall trend of the AHC in these compounds to the changes in their densities of d-states 
around the Fermi energy upon varying the atomic number. Moreover, we show the generality of 
the phenomenon of anisotropic anomalous Hall effect by demonstrating its occurrence within the 
three-band tight-binding model. 
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I. INTRODUCTION 

Despite its long history, the anomalous Hall effect 
(AHE) in ferromagnets, discovered in 1881^ is still not 
fully understood from the theoretical point of view £ Nev- 
ertheless, due to possible vast applications in spintronic 
devices, the AHE, and its counterpart in nonmagnetic 
materials — the spin Hall effect (SHE)^ have drawn 
quite intensive attention in the recent years. The under- 
lying topological nature of the intrinsic AHE and SHE, 
relating these phenomena to some fundamental physical 
effects, makes them even more relevant and interesting. 
Spin-orbit coupling (SOC) plays a crucial role in both 
AHE and SHE, as proposed in the first microscopic the- 
ory of the AHE by Karplus and Luttinger^ It can be 
demonstrated that SOC in perfect crystals gives rise to 
a transverse anomalous velocity of electrons propagating 
along the direction of the external electric field — thus 
leading to the anomalous Hall current. This mechanism 
is nowadays referred to as the intrinsic contribution. 

The intrinsic AHC considered in this work can be ob- 
tained via the linear response Kubo formula for the off- 
diagonal components of the conductivity tensor a: 
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which relates the conductivity tensor to the Brillouin 
zone (BZ) integral of the k-dependent Berry curvature 
tensor f2. In the latter expression and ^ are re- 
spectively the occupied (o) and empty (e) one-electron 
spinor Bloch eigenstates of the crystal, e„k and e m ic are 
their eigenenergies, and and Vj are the Cartesian com- 
ponent of the velocity operator v. As a second-rank an- 
tisymmetric tensor, the AHC tensor can be also seen as 



the anomalous Hall conductivity vector, er, the compo- 
nents of which are related to the components of the AHC 
as <7j = ■^CijkO'jk, where e^k is the Levi-Civita tensor. 

For materials with impurities or disorder, extrinsic 
contributions to the AHE also exist, which can be de- 
scribed within density functional theory££ Nevertheless, 
the intrinsic contribution is often dominating in itiner- 
ant ferromagnets with moderate resistivity^ Since the 
intrinsic anomalous Hall conductivity (AHC) is deter- 
mined by the electronic structure of a pristine crystal 
(Eq. [TJ, which can be accurately calculated using mod- 
ern first principles methods, a comparison between ex- 
periments and first principles calculations serves as the 
first necessary step to deeper understanding of the intrin- 
sic AHE in real materials. Several investigations using 
the first principles methods have been done, for instance, 
in SrRu0 3) Z£ Fe£UL Mn 5 Ge 3 ^- CuCrzSe^Br^ Ni^ 
Co. 13 ' 14 For those materials, the calculated intrinsic AHC 
agrees well with the experimental values, except for the 
case of fee NijiS which might be due to its complicated 
electronic structured 

One of the recently emerging topics in the field 
of the transverse magneto-transport phenomena is the 
anisotropic nature of the off-diagonal part of the conduc- 
tivity tensori 14 ' 16 ' 17 In the case of the AHE, the presence 
of the magnetization M in a ferromagnet leads to a strong 
dependence of the components of the conductivity tensor 
on the magnetization direction in the sample. Although 
experimentally, anisotropic AHE has been observed in 
many materials, e.g. bee Fe^ fee Ni ) 19 ' 20 hep Gdf2i 
as well as FeCr 2 S 4) ^ Yb^MnSbn^ YaFeir-^Co^ 
and i? 2 Fei7 (R = Y, Tb, Gd)d only two studies of 
the anisotropy of the AHE from first principles have 
been performed so far. Roman et alJ^ considered the 
anisotropic AHE in uniaxial hep Co, calculated the ratio 
of the AHCs for the out-of-plane and in-plane magne- 
tization, er z and a x , respectively, and found it to be as 
large as four, which is close to the experimentally ob- 
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served ratioJ^ Moreover, they performed a directional 
averaging of the anisotropic AHC and compared the ob- 
tained conductivity to the experimental value measured 
in polycrystalline hep Co samples^ finding an excellent 
agreement^ In another work, Zhang and co-authorsii 
considered the anisotropic AHE in uniaxial Llo FePt al- 
loy. They also found a large anisotropy of the AHC in this 
compound, and were able to attribute it to the spin- non- 
conserving part of the spin-orbit interaction, prominent 
in this material with strong SOC. 

In this work, we undertake a detailed first principles 
analysis of the anisotropic intrinsic AHE in the group of 
Llo-ordered 3c?Pt (3d = Fc, Co, Ni) alloys. These ma- 
terials are currently under investigation with respect to 
possible spintronic applications due to their large uni- 
axial magnetic anisotropy energies and high Kerr rota- 
tion, making them possible candidates for ultrahigh den- 
sity magnetic and magneto-optical recording media^i 
Recently, the AHE in FePt was used for injection of a 
spin-polarized current for consequent detection of direct 
and inverse spin Hall effect in Au.— In a combined ex- 
perimental and theoretical study it was shown that the 
intrinsic contribution to the anomalous Hall signal dom- 
inates in samples of FePt with finite structural disorder. 
All this motivated our study of the anisotropy of the in- 
trinsic AHE in uniaxial 3dPt alloys. 

In general, we find very large anisotropy of the AHE in 
these compounds, which changes its magnitude and sign 
as a function of the band filling of the 3c? transition- metal. 
In particular, we observe that for the 3c?Pt alloys with 
high concentration of Co atoms the a z and a x conductiv- 
ities differ in sign, which leads to the phenomenon of the 
anti- ordinary AHE, in which at a certain "magic" angle 
of the magnetization the Hall current J becomes parallel 
to M. Moreover, for (Feo.iCo . 9 )Pt and (Co .85Nio.i 5 )Pt 
alloys we predict the occurrence of the colossal anisotropy 
of the AHE, that is, an order of magnitude reduction in 
the value of o~ x as compared to a z , or visa versa. By an- 
alyzing the spin-resolved AHC in these alloys, we try to 
relate the general trend of decreasing AHC in these com- 
pounds when going from FePt to NiPt to the changes in 
their densities of states around the Fermi energy. More- 
over, we demonstrate the occurrence of the anisotropic 
AHE within the "uniaxial" minimal three-band tight- 
binding model, which underlines the generality of this 
phenomenon and hints at its occurrence in a wide range 
of materials. 

The structure of the paper is as follows. In Section II 
we describe the method and details of our first principles 
calculations. In Section III we introduce a minimal three- 
band tig model, necessary to predict the appearance of 
the anisotropy of the AHE in a three-dimensional crystal, 
and investigate the AHC within this model as a function 
of the band filling. In Section IV we present the results of 
our ab initio calculations of the AHE in the family of or- 
dered FePt, CoPt and NiPt alloys. We demonstrate that 
the AHE in these alloys in strongly anisotropic and dis- 
plays a number of interesting phenomena in the region 



[ink A [ooi] 























fpt T 

Fe 





[110] 



FIG. 1: (color online) Crystal structure of Llo FePt alloy. 
Small (red) spheres stand for the Fe atoms, while large (blue) 
spheres mark the Pt ions. The primitive unit cell used in 
the calculations is enclosed with thicker lines. In the text, z 
stands for the [001] axis, while x stands for the [110] direction 
in the crystal. 

where it changes sign. We end the paper with conclu- 
sions. 



II. COMPUTATIONAL DETAILS 

We performed our density-functional theory (DFT) 
calculations of Llo ordered 3c?Pt (3c? = Fe, Co, Ni) alloys 
using the full-potential linearized augmented plane- wave 
(FLAPW) method as implemented in the Julich DFT 
code FLEUR^S The generalized gradient approximation 
(GGA)2I for the exchange-correlation potential was used. 
The self-consistent calculations with SOC were done in 
second variation with k max of 4.0 a.u. -1 and 16000 k- 
points in the full Brillouin zone (BZ). The muffin-tin ra- 
dius of 2.45 a.u. was used for all atoms. Six local orbitals 
for the 4p-states of Pt atoms were used to take care of the 
core charge of Pt correctly. In all our calculations, a tp2 
geometry with two atoms in the Llo phase was used for 
all alloys, with experimental lattice constants (Fig. 1)^ 
For intermediate alloys, for instance, (Fe .5Co .5)Pt, the 
virtual crystal approximation (VCA) was applied on the 
3c? atomic sites, where the composition-averaged core po- 
tential is used instead of that of pure 3d elements, to- 
gether with corresponding number of valence electrons, 
and interpolated lattice constants from the neighboring 
compounds. 

In this work, we applied the Wannier interpolation 
technique to calculate the AHC accurately. We fol- 
lowed the method introduced in Refs. [33[ and [34[ 
to construct the maximally-localized Wannier functions 
(MLWFs) from the FLAPW wave functions, in which 
the unitary transformations are constructed to minimize 
the spread of the Wannier functions. Using the self- 
consistent charge density with SOC included, 36 MLWFs 
corresponding to s,p, c?-orbitals of 3c? and Pt atoms for 
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both spins were generated on a 10 x 10 x 10 fc-mesh, 
using the WANNIER90 code^ The AHC was then calcu- 
lated by applying the Wannier interpolation scheme of 
Wang et ali^ for evaluating the Berry curvature on a 
208 x 208 x 208 uniform fc-mesh. For fc-points at which 
the Berry curvature exceeded 30 A 2 , an adaptive refined 
fc-mesh of 5 x 5 x 5 was used. 



III. AHE ANISOTROPY: GENERALITIES 

In terms of the AHC vector the linear response expres- 
sion for the anomalous Hall current J can be rewritten 
as 

J(M) = <x(M) x E, (2) 

where E is the electric field. In a ferromagnet with uni- 
form magnetization M, er and J depend on the magneti- 
zation direction in the crystal. This magnetocrystalline 
anisotropy of the AHC manifests itself in the changes in 
the direction and magnitude of er upon changing the di- 
rection of the magnetization away from a certain (easy) 
axis. In general, while J is always perpendicular to the 
electric field E, it is not necessarily orthogonal to M, as 
er and M may not be parallel. In single crystals er and M 
are perfectly collinear only when M points along certain 
high symmetry directions. For an arbitrary orientation 
of M there is generally a misalignment between them, 
which is the signature of the anisotropic AHE. Another 
manifestation of the AHE anisotropy is the dependence 
of the absolute value of J and er on the direction of M. 
While in cubic crystals the AHC anisotropy appears only 
at the 3rd order with respect to the directional cosines 
of the magnetization, in uniaxial crystals the linear term 
can dominate^ leading to large AHE anisotropics, ob- 
served experimentally-^ and explained theoretically^ 

The microscopic origin of the anisotropic AHE can 
be easily understood by inspecting Eq. (TTJ). Consider 
a tetragonal crystal structure, as depicted for Ll FePt 
alloy in Fig. 1. Suppose that we start with the magne- 
tization M pointing along the [001] axis (z-axis in the 
following). In this case the v x and v y components of the 
velocity operator have to be inserted at the place of Vi 
and Vj operators in Eq. (jTJ) in order to obtain the a xy 
component of the conductivity tensor, or, cquivalently, 



the a z component of the conductivity vector, taking into 
account that for such a high symmetry direction of the 
magnetization er is aligned with M along the z-axis. Ro- 
tating now M away from the [001] axis modifies (i) the 
wave functions ipnk and ipmk, (ii) occupation of the states 
and (iii) the eigenenergies of the states e„k and e m k due 
to the presence of the spin-orbit interaction. Thus, all 
components of the conductivity tensor have to be recal- 
culated for a general direction of M. In this work, the 
magnetization is confined in the high-symmetry (110)- 
plane, and the resulting conductivity vector also lies in 
the same plane due to the antisymmetric nature of the 
anomalous Hall conductivity with respect to the inversion 
of the magnetization direction. For a general magnetiza- 
tion direction M in this plane, the AHC vector er can be 
decomposed as follows: 

er = (T||M + (Tj_n, (3) 

where M and n are the unit vectors along the magnetiza- 
tion direction and orthogonal to it within the (llO)-planc, 
respectively. The ratio of a||(M) and ctj_(M) gives an es- 
timate of how strongly the AHC vector deviates from the 
direction of M. Upon further rotation the magnetization 
hits the [110] direction in the crystal (x-axis in the fol- 
lowing), and the orthogonal component of the AHC, a±, 
is zero, while er is collinear with the magnetization again. 
In this case v y and v z enter Eq. ([T]), and the magnitude 
of the AHC is given by a x . 



IV. ANISOTROPIC AHE WITHIN THE 

f 29 -MODEL 

In this section we demonstrate the appearance of 
anisotropic AHE within a simple tight-binding model, 
namely, three-band t2 g model for d yz , d zx , and d xy spin- 
up orbitals on a cubic lattice. We consider only the hop- 
pings up to the nearest neighbors, ti, and to the next 
nearest neighbors, £ 2 . The Hamiltonian of the model in 
fc-spacc reads: 

H(k) = H (k)+H so (M), (4) 
where the Hamiltonian without SOC is given by: 



- 2t i (cos k y + A cos fc ; 
Ho(k) = J 4t 2 sin k x sin k y 
4i 2 sin k x sin k z 



, z ) 4£ 2 sin k x sin k y 4i 2 sin k x sin k z 

— 2^ (cos k x + A cos k z ) 4t 2 sin k y sin k z 
4i 2 sin k y sin k z — 2ti(cos k x + cos k y ) 

I 



(5) 



in which we introduced an anisotropy parameter A. The for A ^ 1 the nearest neighbor hopping in the (x, ?/)-plane 
role of this parameter is to make the system uniaxial, i.e., is different from that along the z-axis. In a real cubic 



crystal introducing such a uniaxiality could correspond 
to e.g. changing the interlayer distance along the z-axis 
via application of stress. 

The fc-independent SOC part of the Hamiltonian de- 
pends on the magnetization direction M in the crystal. 
Within our model for M along the z-axis the S yz - and 
dj^-orbitals arc coupled due to SOC, and iJso _ma t r i x 
reads: 



flso(M||z) 




(6) 



while for M along the x-axis d\ x - and d| y -orbitals are 
coupled instead: 













flso(M||a;) = £ I i 
-i 



(7) 



The strength of the spin-orbit interaction is constant in 
both cases and is given by parameter of the model £. The 
band structure of the model obtained by diagonalizing 
Hamiltonian ^ is plotted in Fig. [5] for M||x and A = 1. 

The AHC calculated according to Eq. (jTJ) for the t% g - 
modcl is shown in Fig. [3] as a function of band-filling for 
M||z and M||x, both with A = 1 and A = 0.9. Wc 
assumed for our calculations a lattice constant of 1 A, 
^2 / ii = 0.1 and = 0.02. Obviously, as expected 
from symmetry, the AHC docs not depend on whether 
the magnetization points along the x- or z-axis when 
A = 1, while it displays a strong dependence on the 
electron occupation n, or, equivalently, on the position 
of the Fermi level E F which corresponds to this occu- 
pation. Such sensitive dependence of intrinsic AHC on 
the details of the electronic structure, which stems from 
a very irregular behavior of the Berry curvature in the 
Brillouin zone, is rather well-known.— From our calcula- 
tions it can be seen that when the Fermi energy is posi- 
tioned in the vicinity of the band edges with high density 
of electronic states (DOS) — which in our model corre- 
sponds to the case of nearly filled, half-filled and com- 
pletely filled bands — the Ep position has a strong effect 
on the AHC (Fig. [3]). The reason behind such a sensitive 
dependence lies in the presence of flat degenerate bands 
around the Fermi energy, which provide wide regions in 
fc-space where the occupied and unoccupied states are 
separated by a small energy. We speculate that such a 
situation can also be related to the anomaly of the den- 
sity of states near the band edges and associated Lifshitz 
transitions^ We expect that at such transitions, anoma- 
lies in the AHC would lead to e.g. singular behavior of 
the anomalous thermopower. 

In case of the "uniaxial" ^g-model (A = 0.9) the elec- 
tronic structure of the crystal with M||x and M||z is not 
the same anymore, and a strong anisotropy of the AHC 
can be seen in Fig. [31 leading even to a difference in sign, 
considered in detail in the following section. In analogy to 
the case of the "isotropic" crystal, the largest difference 
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FIG. 2: Electronic bands of the t2 g -model for M||a:, with 
t2/t\ = 0.1, = 0.02, A — 1.0, and the lattice constant 
of 1 A. The dashed horizontal lines mark the position of the 
Fermi level for the electronic occupation given by the number 
on the right. Labels mark the orbital character of the bands. 



between the AHC for the two magnetization directions 
can be observed for n » 0, n « 1.5 and n w 3, although 
for exact half-occupation (n ~ 1.5), this difference van- 
ishes. 

In order to see the origin of the anisotropic AHC in the 
uniaxial <2g-crystal clearly, we refer to the band structure 
of the system, Fig. [21 for M||x and A = 1. Let us consider 
the case of n = 1.2 and transitions between the occu- 
pied and unoccupied states in the energy region marked 
with red solid circle along the TA-path. For M||x the 
SOC leads to the mixing of d xy - and <i z:E -orbitals, and 
the resulting small energy splitting between the two cor- 
responding bands can be clearly seen. The non- vanishing 
matrix element of the SOC between the two latter or- 
bitals leads to a finite contribution to the Berry curvature 
and the AHC, associated with the electronic transitions 
across E F 2 . On the other hand, the d zx - and d^-orbitals 
are not coupled by SOC for this magnetization direction 
(which can be also seen from an exact degeneracy of cor- 
responding bands along TX'), and the contribution from 
the symmetry-equivalent without SOC part of the band 
structure along FX' to the Berry curvature, marked with 
dashed blue circle, is exactly zero. The situation is re- 
versed for M||z in the isotropic crystal, and the contri- 
bution to the AHC from the states in the dashed blue 
circle is exactly the same as that from the corresponding 
region along the TA-path for M||x, while the latter gives 
no contribution for M||z. Overall, when only the encir- 
cled regions and their symmetric "clones" are considered, 
the resulting AHC does not depend on the magnetization 
direction and there is no anisotropy of it. Introducing 
now anisotropy in the system by setting A to 0.9 in the 
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FIG. 3: Anomalous Hall conductivity as a function of the 
band filling within the three-band t2 S model with the param- 
eters given in the caption to Fig. [5] Open (filled) circles mark 
the case of M||x(z) in the "isotropic" crystal with A = 1 
in Eq. Q. Note that in this case the AHC curve for M||z is 
shifted by 5 S/cm with respect to the AHC with M||x in order 
to see the degeneracy between the two curves clearly. Filled 
squares (diamonds) stand for the AHC in the "anisotropic" 
crystal with A = 0.9 in Eq. © for M||x(.z). 



i2 g -model leads to the fact that the electronic structure 
along the FX- and TAT'-paths is not the same anymore 
and thus, the contributions to the AHC from the full red 
circle for M||a; and from the dashed blue circle for M||x 
are different. In turn, this leads to the AHE anisotropy. 
This line of thinking is clearly valid also for explaining 
the anisotropy of the total AHC, which is obtained by a 
summation over all such encircled regions in energy and 
fc-space, contributing to the AHC. 



V. ANISOTROPIC AHC IN 3dPt ALLOYS 

The results of our calculations for the intrinsic AHC in 
Llo-ordered FePt, CoPt and NiPt and their intermediate 
alloys are presented in Tableland in Fig.@]for [001] (tr z ) 
and [110] (a x ) directions of the magnetization M in the 
crystal. In general, the observed behavior of a z and a x as 
a function of the electron occupation of the 3d transition- 
metal is similar: starting from the FePt alloy with posi- 
tive values of the AHC for both magnetization directions 
of around 600 S/cm, the sign of o~ z and a x changes in 
the vicinity of the CoPt alloy, and the AHC values are 
very large and negative for NiPt, reaching as much as 
— 1200 S/cm. On average, we can see that the decrease 
of both conductivities with increasing electron occupa- 
tion is rather linear. It is worth mentioning that such 
decreasing trend of the AHC is somewhat reminiscent of 
the trend among the pure bec Fe, hep Co and fee Ni, for 
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TABLE I: Values of the AHC in Llo FePt, CoPt and NiPt 
with the magnetization along [001] (o z ) and [110] (o x )- For 
each orientation, o^ (o^) is calculated by keeping only the 
first (second) term in the spin-orbit Hamiltonian (Eq. [SJ, 
while both terms are kept when calculating <r tot . Ao tot is 
defined as o z —o x - Acr 1 "" 1 " 1 "' is defined as the difference between 
the spin-conserving (spin-flip) parts of o z and o x . All values 
are in S/cm. 



which the calculated intrinsic values of the u z AHC con- 
stitute approximately 750 S/cm^ii^ 480 S/cm / 3 ' 14 and 
—2200 S/cm^ respectively. We will come back to this 
point at the end of this section. 

It is clear from Table Q] and Fig. 2] that for almost all 
considered alloys the anisotropy of the AHC reaches very 
large values. This is expected, since in uniaxial crystals 
the AHC anisotropy appears already in the first order 
with respect to the directional cosines of the magnetiza- 
tion (see discussion in the previous section)^ In FePt, 
the difference between a z and a x , Aa tot (filled circles in 
Fig. is as large as a x itself and constitutes around 
400 S/cm, c.f. Table Qp! In CoPt, on the other hand, 
the absolute value of Acr tot is twice larger than the ab- 
solute value of the AHC for any of the two magnetiza- 
tion directions. The AHC anisotropy reaches as much as 
—500 S/cm in the vicinity of FeCoPt and CoNiPt alloys, 
and in general, the behavior of Acr tot is neither smooth 
nor monotonous, displays several mimina and maxima as 
a function of the electron occupation, and even changes 
its sign for Fe^Coi-zPt alloy with x « 0.75. On the other 
hand, the anisotropy of the AHE in 3dPt alloys, when 
the magnetization is rotated in the (OOl)-plane, is much 
smaller than the "out-of-planc— in-plane" anisotropy dis- 
cussed previously, which can be easily understood taking 
into consideration the higher symmetry of the former sit- 
uation. E.g., the difference of the AHCs for M along [110] 
and [100] reaches at most 80 S/cm in NiPt alloy, being 
as small as —47 S/cm in CoPt and —16 S/cm in FePt. 

With grey shaded area in Fig. 2] we highlight the 
region around the CoPt alloy, where both a z and o~ x 
change their sign. This sign change leads to the oc- 
currence of two interesting phenomena with respect to 
the anisotropic AHE. The first one, which we name the 
colossal anisotropy of the AHE, according to our calcu- 
lations, occurs for Fe^Coi-^Pt alloy with x ~ 0.1 and 
for Co^Nii-^Pt alloy with x w 0.85. For these two com- 
pounds one of the conductivities, a z for Feo.1Coo.9Pt and 
a x for Coo.85Nio.15Pt, turns to zero, which marks the 
complete disappearance of the intrinsic anomalous Hall 
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FIG. 4: Anomalous Hall conductivity of 3dPt alloys for M 
along [001] (<Tz, open circles) and [110] {a x , open squares), 
and anisotropy (A<j iot = a z — a x , filled circles) , with respect 
to the band filling. 

current J for one of the magnetization directions in the 
crystal. We introduce the term colossal anisotropy in 
analogy to the situation which was predicted to occur 
in one-dimensional Pt wires, for which upon changing 
the magnetization direction M the value of the mag- 
netization |M| itself can be quenched completely— In 
terms of the longitudinal transport within the setup 
of e.g. anisotropic magnctorcsistancc (AMR) experiment, 
the occurrence of the colossal anisotropy of the diagonal 
conductivity would results in a metal-insulator transition 
in the crystal — in case of the colossal AHE anisotropy 
observed in 3dPt alloys all compounds remain metallic for 
all magnetization directions, however, and retain their 
complicated electronic structure around the Fermi en- 
ergy 

For Coo.8Nio.2Pt alloy in Fig. [5ja) we plot the depen- 
dence of the cry and <r± components of the AHC on the 
angle 9 of the magnetization with the z axis when it is 
rotated away from the [001] direction within the (110) 
plane towards the [110] direction. At 9 = the a± com- 
ponent is zero and the AHC vector with magnitude of 
340 S/cm is antiparallel to the z axis, along which the 
magnetization is aligned, c.f. Fig. [SJc). Upon increasing 
9 we observe the increase in cry and decrease in <jj_, with 
both components becoming equal at the angle 9 ps 55°. 
At this angle the magnitude of the AHC is reduced signif- 
icantly to 210 S/cm, while its deviation from the z-axis is 
only about 10°. Thus, in this range of 9, the rotation of 
the magnetization results mainly in quenching the mag- 
nitude of the anomalous Hall current, while its direction 
basically remains "stuck" to the [110] axis. Upon fur- 
ther rotation of the magnetization both components of 
the AHC vector rapidly approach zero, the AHC vector 
quickly rotates towards the — x axis, and when M hits 
the [110] direction, the AHC with a tiny magnitude of 
25 S/cm is again antiparallel to the magnetization. 




FIG. 5: (a) Colossal anisotropy of the AHC in Coo.8Nio.2Pt 
alloy. Red circles (blue squares) denote the erii (tJ±) compo- 
nent of the AHC as a function of the angle 9 of the mag- 
netization M with [001]-axis upon rotating it into the [110] 
direction, (b) Anti-ordinary Hall effect in CoPt. Red circles 
(blue squares) denote the cry (cr±) component of AHC, as a 
function of the angle 6 of the magnetization M with [001]- 
axis upon rotating it into the [110] direction, (c)-(e) depict 
the relative orientation of the Hall current J, AHC a and 
magnetization M in the situation of the anti-ordinary AHE. 
In (c)-(e) the magnetization is confined to the (llO)-plane. 



For CoPt alloy the situation, depicted in Fig.[5jb-e), is 
completely different. Similarly to the previously consid- 
ered case, at 9 = 0° the AHC vector is antiparallel to M, 
and its magnitude constitutes 120 S/cm, Fig. [Sic). Upon 
increasing 9 up to as much as 45° the conductivity vec- 
tor resides basically in the close vicinity of the [001] axis, 
while its magnitude increases. For example at 6 = 45°, 
(7|| « cr j_ and the value of total a is roughly 170 S/cm, 
Fig. [5jd). With further increasing 9 the magnitude of 
the AHC is increasing even further, while the AHC vec- 
tor starts its way towards the [110]-direction. The in- 
crease of I cr I is mainly due to the a± component in this 
regime, while at the same time |<T||| is becoming smaller, 
and eventually changes its sign. Finally, at 9 = 90°, the 
AHC vector is aligned together with M along the x axis, 
and its magnitude is 110 S/cm. 

Remarkably, cry turns to zero at 9q — 70°, which man- 
ifests the occurrence of the anti- ordinary Hall effect in 
the crystal of CoPt, see Fig.[SJc). At this "magic" angle, 
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the magnitude of the anomalous Hall current J is almost 
twice larger than it is for M||z, however, due to non- 
vanishing o~± component of the AHC vector, J is aligned 
along the direction of the magnetization. By analyzing 
Figs. EJb-e) we observe that the rotational sense of the 
anomalous Hall current is opposite to that observed in 
the ordinary Hall effect (OHE) of free electron gas. For 
OHE, Lorentz forces ~ [H x v] are acting on electrons 
with velocity v in the presence of magnetic field H. The 
resulted ordinary Hall current of free electrons is always 
perpendicular to H irrespective of its direction, opposite 
to the situation of the anti-ordinary anomalous Hall ef- 
fect, observed in CoPt. Here, turning the magnetization 
clockwise in the (llO)-planc results in an anti-clockwisc 
rotation of J, with its value staying rather large all the 
time. The anti-ordinary spin Hall effect has been also 
recently predicted to occur in transition metals.— 

In the region of 3dPt alloys in the vicinity of Llo CoPt 
the anisotropy of the AHE manifests itself in crucial ways 
suggesting new functionalities of the AHE-based devices. 
In this region, large changes in the magnitude of the 
anomalous Hall current as well as relative orientation of 
the Hall current with respect to the magnetization can 
be easily achieved by simple reorientation of the sam- 
ple's magnetization. While the former could be used in 
order to e.g. tune the relative magnitudes of the extrin- 
sic and intrinsic anomalous Hall signal^^ among most 
straightforward applications of the latter could be a re- 
alization of the planar Hall effect (PHE)^ which is re- 
lated to the Hall effect in ferromagnetic materials ob- 
served in a two-dimensional geometry with electric field, 
magnetization and the Hall current sharing same sample 
plane. So far, it is believed that in most of the cases 
the PHE originates from anisotropic magnetoresistance 
in metallic ferromagncts, although the PHE mechanism 
stemming from the anomalous Hall effect due to non- 
collinearity of the magnetization in semiconductor-based 
materials has been also suggested.— Within the scope 
of the anti-ordinary Hall effect, described in this work, 
it would be possible to observe the PHE coming solely 
from the anisotropic nature of the collinear ferromagnetic 
materials. 

We would like to underline, that despite the crude- 
ness of the VCA approximation for description of the 
electronic structure of complex alloys, the results of our 
work still hold, although the exact width of the region 
where the colossal anisotropy and anti-ordinary nature 
of the intrinsic anomalous Hall effect can be observed, 
might be different when more appropriate approxima- 
tions, such as coherent potential approximation (CPA)fS 
are used to treat the substitutional alloys Fe^Coi-a-Pt 
and Co^Nii-^Pt. The main reason behind this is that 
for "pure" ferromagnets FePt, CoPt and NiPt our re- 
sults are exact in the sense that no approximations of 
disorder need to be made and the AHE consists only of 
the intrinsic contribution, while the precise value of the 
intrinsic AHC will still depend on the chosen parame- 
ters and formulations of the DFT calculations such as 



exchange-correlation functionals, treatment of SOC, va- 
lidity of the single-particle picture, particular choice of 
the basis set, etc. The values of the intrinsic AHC at 
the ends of the considered family of alloys, namely, FePt 
and NiPt, are large in their magnitude but differ in their 
sign. This means that upon varying the concentration x 
in Fe^Coi-^Pt and Co^Nii-^Pt alloys, the region where 
the AHC changes sign must exist, irrespective of the ap- 
proximations made. At the end, it is the change of sign of 
the AHC for the CoPt alloy which leads to the occurrence 
of the colossal anisotropy and anti-ordinary anomalous 
Hall effect in its vicinity according to our calculations. 

At the end of this section we will try to relate the 
mentioned above change of sign in the values of the AHC 
when going from FePt to NiPt, to the changes in the elec- 
tronic structure of these materials. For this purpose, we 
first of all decompose the atomic spin-orbit Hamiltonian 
in the well-known way>^£ 

£L ■ S = £L fi S fl + C (i+ST + L7St) /2, (8) 

where £ is the spin-orbit coupling strength, h is the 
spin magnetization direction (which is taken as the spin- 
quantization axis) , L and S are the total orbital and spin 
angular momentum operators, Lf, = L • n, and Lt and 
Lr are the corresponding raising and lowering operators 
(analogously for spin). We shall refer to the first and 
second terms in Eq. (jSJ) as the spin-conserving and spin- 
flip parts of the SOC. This terminology refers to the ef- 
fect of acting with each of them on an cigenstate of Sa. 
Accordingly, we define and as the AHC calcu- 
lated from Eq. (1) after selectively removing the second 
or the first term on the right-hand side of Eq. (2). This 
is not an exact decomposition, but inspection of Table U 
shows that it is approximately valid for both magneti- 
zation directions in FePt and CoPt, and NiPt with M 
along x, i.e. a tot w er^ in these cases, while even for 
NiPt with M||z a large discrepancy between the er tot of 
— 1165 S/cm and the sum + a'®' ps —2000 S/cm does 
not change the general line of argument we are to present 
below. 
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TABLE II: Spin-resolved contributions to the spin-conserving 
AHC in 3dPt alloys. SOC 3d (Pt, 3d+Pt) stands for the val- 
ues obtained with SOC on only 3d (Pt, both 3d and Pt) site(s) 
included in the calculations, "f (4.) denotes the contribution 
from the majority (minority) spin channel. All values are in 
S/cm. 

By analyzing Table U we observe that the spin- flip con- 



ductivity in 3<iPt alloys provides a significant contribu- 
tion to the total AHC, which is particularly striking in 
case of CoPt where is even somewhat larger than the 
spin-conserving part for both magnetization directions. 
And while in NiPt the AHC anisotropy Atr tot is mainly 
given by the anisotropy of its spin-conserving part, Acr^, 
in FePt and CoPt the anisotropy of the AHC is driven 
entirely by the anisotropy of c^-, Ac^-, which exceeds as 
much as 90% of Acr tot in FePt. Such a pronounced role 
of the spin-flip SOC for the anomalous Hall conductiv- 
ity and its anisotropy in ferromagncts containing heavy 
elements, such as Pt, was demonstrated and explained 
by Zhang and co-workers by employing the perturbation 
theory arguments^ In the case of considered here al- 
loys, from Table Hit is however clear that, despite a large 
spin-flip contribution, the overall trend of the total AHC 
between FePt and NiPt can be qualitatively described by 
considering the spin-conserving AHC only, and we dedi- 
cate the rest of the paper to the analysis of <r^ in FePt, 
CoPt and NiPt compounds. 

Firstly, the advantage of considering exclusively the 
spin-conserving SOC is that spin remains a good quan- 
tum number, and the conductivity can be unambigu- 
ously decomposed into spin-up and spin-down parts: 
0-TT = -t_)_ cr 4-_ j n particular, this means that if we assume 
that in the system considered the spin-orbit is given only 
by the spin-conserving part, the corresponding spin Hall 
conductivity ctsh can be obtained as the difference be- 
tween the spin-resolved conductivities: crgfj = <J^ — cp , 3 
which implies that in a non-magnetic material, such 
as e.g. Pt, cr^ = —cr - , while ctsh = — 2cp. Secondly, 
among the cr" and cp* conductivities, the later one is 
much more sensitive to the details of the Fermi sur- 
face, while the electronic transitions contributing to the 
spin-conserving AHC according to Eq. (1) arc distributed 
much broader in energy around Ep^ which makes the 
analysis of the latter easier. 

In order to get additional insight into the structure of 
cr^, we use the atomic decomposition of the AHC for 
each spin channel, considered by Zhang et al.^ based 
on the following atomic decomposition of the spin-orbit 
part of the Hamiltonian: 

tfso=e3dL 3d -S + £p t L pt .S, (9) 

where L 3d ( pt ) is the orbital angular momentum opera- 
tor associated with 3d (Pt) atoms, and ^3d(pt) is t ne 
spin-orbit coupling strength averaged over the valence 
d-orbitals inside 3d (Pt) atom, with the values of 0.54 eV 
for Pt and 0.05—0.07 eV for 3d transition-metal atoms. 
By selectively turning off the spin-orbit coupling inside 3d 
transition-metal atoms (£ 3£ ; = 0) or Pt atoms (£p t = 0) 
we obtain the values of <Tp^ and a^j , respectively. 

The results of our calculations for the spin and atomi- 
cally decomposed in FePt, CoPt and NiPt alloys are 
presented in Table [TTJ Let us take a look at the first two 
columns of the table, where the values of the total cp and 
cp are listed. Firstly, we observe that positive cp and 
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E - E F (eV) 

FIG. 6: Atomically-resolved density of d-states in FePt, CoPt 
and NiPt alloys. Up- and down-arrows stand for spin-up and 
spin-down. 



negative cp arc opposite in their sign for all alloys. Sec- 
ondly, upon going from FePt to NiPt, the spin- up AHC 
increases but retains its order of magnitude, being about 
650 S/cm for FePt and 1300 S/cm for NiPt. On the other 
hand, a very small spin-down AHC of «— 100 S/cm for 
FePt increases by an order of magnitude and reaches as 
much as —2500 S/cm in NiPt. Correspondingly, in FePt 
the positive sign of <r^ is due to the AHC in spin- up 
channel, while in NiPt the negative cr^T is driven by large 
and negative spin-down AHC. 

Consider now the case of FePt. The atomic decom- 
position of the AHC, presented in Table [H] clearly re- 
veals that the large spin-up AHC in this alloy originates 
from the spin-up contribution of Pt atoms, while the Fe 
contribution to cr^ is very small. In the spin-down chan- 
nel, Pt and Fe AHCs, both with the magnitude of about 
200 S/cm, are opposite in sign and suppress each other. 
In CoPt the spin-up Pt and Fe AHCs remain basically the 
same compared to FePt, while the corresponding spin- 
down conductivities significantly increase in their magni- 
tide. This can be related to the increase in both Co and 
Pt density of states of d-electrons around the Fermi en- 
ergy for minority spin which can be clearly seen in Fig. [6j 
as compared to respective DOS of FePt alloy. Such an 
enhancement of the DOS around Ep results in more oc- 
cupied and unoccupied d-states and corresponding tran- 
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sitions accross the Fermi energy, which contribute to the 
AHC according to Eq. (1) — similar to the situation we 
came accross when analyzing the tight-binding model re- 
sults previously. While in the latter case the variations of 
the Fermi energy in the region of increased DOS resulted 
in large changes of the AHE, in the case of a complex 
ferromagnet with many bands at Ep in which the AHE 
is not driven by a single band degeneracy, it seems rea- 
sonable to assume that the increased number of available 
transitions will lead to a larger magnitude of the AHC. 
The increased Co and Pt spin-down AHCs are still oppo- 
site in sign however, which still suppresses the total cr^, 
although its value is also somewhat enhanced compared 
to FePt, which leads to the decrease in overall a = cA+cr^ 
in CoPt, c.f. Table I. 

The reason behind increased d-DOS of Pt atoms at 
the Fermi level in CoPt lies in moving of the spin-down 
Co (i-subband to lower energies with decreasing exchange 
splitting. This leads to a stronger hybridization between 
the Co and Pt d-states, which are situated mainly be- 
low Ep, and increased Pt DOS (Fig. [H]). This is even 
more pronounced in case of NiPt, where the spin-down Ni 
subband lies predominantly below the Fermi energy and 
the hybridization with the Pt d-states is even stronger 
(Fig. [5]) . Correspondingly, as a result of even more en- 
hanced spin-down d-DOS of Ni and Pt atoms around Ep 
in NiPt alloy, the values of tr^j and Op t become very large, 
reaching as much as —1600 S/cm for Ni and —2300 S/cm 
for Pt. On the other hand, the increase in Ni and Pt 
AHC in the majority channel is quite moderate due to 
slightly enhanced DOS, and the total AHC in NiPt be- 
comes large and negative. 

The change of sign of the 3d spin-down subband AHC 
between Fe, Co and Ni in 3dPt alloys can be probably 
related to different orbital character of the d-states at 
the Fermi energy and corresponding matrix elements of 
the SOC in these transition metals, which could in turn 
explain the AHC sign change in elemental Fe, Co and Ni, 
observed experimentally, and reproduced from the first 
principles^ ' 10 ' 13 i 14 Such a change of sign as a function 
of the Fermi level position within the spin subband has 
been also demonstrated from first principles calculations 
of the spin Hall conductivity in Pt^S as well as from 
tight-binding calculations of the spin Hall conductivity 
in 4d and 5d transition-metals4i Assuming that in latter 
cases the spin-flip contribution to the spin Hall effect is 
negligible^ this results in a corresponding sign changes 
of the Hall conductivity, discussed previously. 

From Table II we can see that the AHC originating 
from the Pt atoms is generally larger in magnitude than 
that from 3d transition-metal. This can be explained 
by noticing that the large spin-orbit constant inside Pt 
atoms £pt is by an order of magnitude larger than £3^. 
Moreover, a consistently positive and negative sign of 



large a\ t and (T Pt throughout the 3dPt family can be re- 
lated to a small spin-polarization of the Pt atoms, i.e. un- 
der the condition that this spin-polarization is at all ab- 
sent, Op t should be equal to — Op t , while the difference 
of the both would provide a value of the intrinsic spin 
Hall conductivity in Pt of about 2000 S / cm. As we can 
see from Table II, these arguments indeed explain the 
sign and magnitude of the Pt-originated AHC. Finally, 
we would like to remark that although the atomic de- 
compositon for the spin-resolved conductivities is over- 
all rather reasonable in that « Cp t + a^ d and sa 
<7p t + Cg d , see Table II, such a decomposition works much 
better for the majority channel. We attribute this obser- 
vation to a much stronger hybridization of the 3c? and 
Pt d-states for the minority-spin around the Fermi level, 
which enhances the contributions to the AHC for which 
the presence of SOC on both Pt and 3d transition-metal 
atoms is important. Such contributions arc omitted in 
the atomic decomposition used above. 



VI. CONCLUSIONS 

In conclusion, we investigated from the first principles 
the intrinsic anomalous Hall effect in 3dPt alloys. From 
our calculations it follows that the AHC in this type 
of compounds is strongly anisotropic. We demonstrate 
the generality of such anisotropy in uniaxial ferromag- 
nets by considering a simple three-band tight-binding 
model. In combination with the sign change of the con- 
ductivity upon going from FePt to NiPt the pronounced 
AHC anisotropy leads to the occurrence of the colossal 
anisotropic AHE and anti-ordinary AHE in the vicinity of 
the CoPt alloy. While in the case of colossal anisotropic 
AHE the anomalous Hall current completely vanishes for 
one of the magnetization directions in the crystal, within 
the scope of the anti-ordinary AHE the rotational sense of 
the Hall current is opposite to that of the magnetization, 
and a complete collinearity of the two can be achieved 
for a certain "magic" angle of the magnetization in the 
crystal. We relate the general trend of the AHC in these 
alloys to the changes in their electronic structure in the 
vicinity of the Fermi level, and discuss possible applica- 
tions of the anisotropic AHE in these compounds. 
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